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Electron-Hole Asymmetry in Superconductivity of Pnictides Originated from the 
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We have performed a systematic photoemission study of the chemical potential shift as a function 
of carrier doping in a pnictide system based on BaFe2As2. The experimentally determined chemical 
potential shift is consistent with the prediction of a rigid band shift picture by the renormalized 
first-principle band calculations. This leads to an electron- hole asymmetry (EHA) in the Fermi 
surface (FS) nesting condition due to different effective masses for different FS sheets, which can be 
calculated from the Lindhard function of susceptibility. This built-in EHA from the band structure, 
which matches well with observed asymmetric superconducting domes in the phase diagram, strongly 
supports FS near-nesting driven superconductivity in the iron pnictides. 

PACS numbers: 74.70.Pq, 71.30.+h, 79.60.-i 



Even though the detail of the pairing mechanism in 
the recently discovered iron-based superconductors is still 
under intense debate, several theoretical investigations 
[I-S and experimental observations [6-10 strongly sug- 
gest the importance of inelastic inter-band scattering 
between hole and electron Fermi surface (FS) pockets 
connected via the antiferromagnetic (AF) wave vector. 
Within this framework, the pairing strength depends on 
near- or quasi- nesting, here defined as a large enhance- 
ment of the spin susceptibility at a well defined wave vec- 
tor |5J . The near- nesting conditions depend on the shape 
and size of the various FS pockets, which are tuned by 
the position of the chemical potential. The evolution of 
the chemical potential with carrier concentration is thus a 
key issue to understand the evolution of FS near-nesting 
and superconductivity in these materials. 

The 122-structural phase of BaFe2As2 is particularly 
suitable for a systematic study of the chemical potential 
shift since it can be doped either by electrons or holes fol- 
lowing Fe^+^ Co^+ or Ba^+^ K+ partial substitutions, 
respectively. Interestingly, the electron- and hole-doped 
sides of the phase diagram show some noticeable differ- 
ences. For example, while the maximum Tc value for 
the hole-doped side reaches 37 K at ambient pressure, it 
tops around 25 K for the electron-doped systems. Simi- 
larly, the superconducting dome extends to much higher 
doping in the hole-doped case, with an optimal concen- 
tration of around 0.2 hole/Fe against 0.08 electron/Fe for 
the electron-doped side. 

Angle-resolved photoemission spectroscopy (ARPES) 
is a powerful tool to access directly the electronic struc- 
ture with respect to the chemical potential. Our previous 



ARPES studies have already revealed strong variations 
in the pairing strengths associated with the various FS 
sheets in the electron-doped compounds [9 as compared 
to the hole-doped ones [6, 11 , as well as the deterioration 
of the near-nesting conditions in highly overdoped sam- 
ples for which Tc vanishes or is significantly suppressed 
[HI [12] . Although near- nesting was naturally proposed to 
explain these anomalies, this concept has not been linked 
to the origin of the electron-hole asymmetry (EHA) and 
up to date there is still no systematic investigation of 
the impact of the chemical potential shift on the band 
structure throughout the whole phase diagram. 

In this letter, we present a systematic ARPES study 
of the chemical potential as a function of carrier dop- 
ing in the 122-pnictides. With doping, the chemical po- 
tential moves smoothly with respect to the low-energy 
valence states, in agreement with our local density ap- 
proximation (LDA) calculations. However, we observed 
anomalously larger (smaller) core level shift than the va- 
lence band shift on the hole (electron) -doped side for the 
relatively undisturbed As ?>d core levels, possibly due to 
the screening effect which increases (decreases) the core 
level shifts upon hole (electron) doping. Based on a rigid 
band shift approximation justified by our experimental 
results, we computed the doping dependence of the Lind- 
hard spin susceptibility at the AF wave vector, and found 
that the Lindhard function itself is asymmetric as a func- 
tion of doping, in a similar fashion as the asymmetry 
between the hole and electron superconducting domes. 
This strongly supports FS-near-nesting-enhanced super- 
conductivity in the pnictides. 

The high-quality single crystals of the 122 series used 
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(Tc = 37 K), and Bao.aoKo.roFesAss (T^ = 22 K). For 
simplicity, here after we call them C00.30, Coo.ie, Coq.os, 
Bal22, Ko.25, K0.40 and Kq.to, respectively. From low 
to high binding energies, we observed the Fe Zd (around 
the Fermi level), Ba 5p (-14.5 eV), K 3p (-18 eV), As 
U (- 40.4 and 41.3 eV) and Ba M (- 89.5 and 92 eV) 
states, respectively. In particular, the As Zd peaks are 
very strong in all compounds regardless of Co and K 
doping. Based on a previous photoemission study ^ 
that indicates that the As atoms in BaFe2As2 are not 
perturbed significantly at the cleaved surface, we used the 
As Zd core levels to investigate the doping dependence of 
the chemical potential. In Fig. 1(b), we show a zoom of 
the As Zd core levels of all compounds. The position of 
the peaks moves towards the lower binding energies as K 
concentration increases. In contrast, the peak positions 
are almost unaffected by Co-doping. We plot in Fig. 1(c) 
the shift of the As 3^3/2 and As 3^5/2 levels as a function 
of carrier density, which is half of the value x for both K 
and Co dopings. The blue dashed line in Fig. 1(c) gives 
the average of the As Zd peaks as a function of doping. 



FIG. 1: (a) Core levels of the Bal22 series recorded with 
a photon energy of 140 eV. The inset shows a zoom of the 
core levels of the K0.70 sample in the 12-21 eV binding energy 
range, (b) Zoom of the As 3d core levels, (c) Doping depen- 
dence of the As 3(^3/2 and As 3(i5/2 core level energies as a 
function of carrier density (half of the x value) . The average 
is represented by the blue dotted line. 



in this study were grown by the fiux method [T3^. Low- 
energy electron diffraction on mirror-like cleaved surfaces 
show a sharp 1x1 pattern in the non-magnetic phase. 
High-resolution (4-20 meV) ARPES measurements of the 
low-energy electronic structure were performed in the 
photoemission laboratory of Tohoku University using a 
microwave-driven Helium source (hv = 21.218 eV) and 
core level studies were done at the Synchrotron Radia- 
tion Center and the Advanced Light Source in USA, as 
well as at the Photon Factory in Japan, using various 
photon energies. Our experiments have been performed 
using high-efficiency VG-Scienta SES-100, SES-2002 and 
R4000 electron analyzers. Samples were cleaved in situ 
and measured at 7-40 K in a vacuum better than 1 x 
10~^^ torr. The samples have been found to be very sta- 
ble and without degradation for the typical measurement 
period of 20 hours. 

Photoemission allows measurement of the core level 
states relative to the chemical potential. It has been used 
widely in the past to study the chemical potential shift in 
high- Tc cuprates [14-16 . Fig. 1(a) shows a comparison 
of the core levels in the 0-110 eV binding energy range 
of 7 samples distributed in the electron-doped and hole- 
doped sides of the phase diagram. These samples are 
BaFei.7oCoo.3oAs2 (T, = K), BaFei.84Coo.i6As2 (T^ = 
20 K), BaFei.92Coo.o8As2 (T, = K), Ba2Fe2As2 (T^ = 
K), Bao.75Ko.25Fe2As2 (T^ = 26 K), Bao.6oKo.4oFe2As2 
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FIG. 2: (a) and (b) ARPES intensity plots of the Coo.ie and 
Ko.40 samples, respectively, along the cut passing through M 
indicated by a dashed line in panel (e). (c) and (d) Corre- 
sponding second derivative intensity plots, (f) Bottom of the 
electron bands versus carrier density. 



An alternative and more direct determination of the 
chemical potential shift is obtained by looking at the 
band dispersion near the Fermi level {Ep)- In Figs. 2(a) 
and (b), we present ARPES intensity plots of the Coo.ie 
and Ko.40 samples along a cut passing through M as in- 
dicated in Fig. 2(e). The corresponding second deriva- 
tive intensity plots are displayed in Figs. 2(c)-(d). The 
blue dashed lines are guides to the eye indicating the 
bottom of the upper electron band (the 7 band as de- 
fined in Ref. [T8|) . The bottom of this electron band at 
the M point moves down from Ep as the signed concen- 



tration decreases (more electrons), which is what we ex- 
pect from simple band filling. In particular, this behavior 
supports the assumption that the Fe -^ Co substitution 
electron-dope the Fe layer, in contrast to a recent density 
functional theory calculation suggesting that Co and Ni 
only act as scattering centers in the Fe planes [19]. It 
is also consistent with the observation of a downshift of 
the F-centered holelike bands in the Co-doped side [8 . 
Fig. 2(f) summarizes our results of the seven differently 
doped samples and gives the position of the bottom of 
the electron band as a function of the carrier density. 
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FIG. 3: Valence band and core level shifts as a function of 
the carrier density, (a) The black line with open circles shows 
the core level shifts while the blue line with closed circles gives 
the valence band shifts measured from the band bottom. The 
red line with open squares gives the difference between the 
core level and chemical potential shifts, (b) The blue dots are 
the valence band shifts shown in panel a, and the red dash 
line is the LDA calculated values of the chemical potential 
divided by a factor of 4. (c) Pictorial representation of the 
explanation of the core level and the chemical potential shifts 
as a function of carrier density. 



At this point, it is instructive to compare LDA calcu- 
lations to the core level shifts and the shift of the bottom 
of the electron band, which corresponds to the chemical 
potential shift in a rigid band picture. The results are 
summarized in Fig. 3. It is clear from Fig. 3(a) that the 
core level shift is not the same as the shift of the valence 
band, and the difference between them are larger on the 
hole-doped side, which will be discussed below. Inter- 
estingly, the theoretically calculated chemical potential 
shift is very much consistent with the observed valence 
band shift when theoretical values are divided by 4 as 
shown in Fig. 3(b), which is understood in terms of the 
band renormalization reported in previous ARPES stud- 
ies H m [18]. This indicates that the shift of the 
valence band corresponds to the chemical potential shift, 
and consequently, the rigid band picture derived from the 
renormalized band structure is valid. 

The core level shift can be understood as follows. The 
core level shift A^ is related to the chemical potential 
shift A/i by the relation: 

AE = -A/i + KAQ + AVm + AEr 

where AQ is the change in valency, K is a constant, 
AVm is a shift due to change in the Madelung potential. 



and AEr is the change in the core-carrier screening [20] . 
Doping is not expected to change the As valency. This 
implies that the term KAQ can be neglected. Therefore, 
the difference between the core level and the chemical 
potential shift represented in Fig. 3(a) by the red line 
is only related to AVm and AEr. It is known that the 
screening term AEr is proportional to the mobile carrier 
concentration, thus one expects that it has the same sign 
on the electron- and hole-doped sides and increases with 
doping. Such doping dependence of the screening term, 
as indicated in Fig. 3(c), will increase (reduce) the core 
level shift caused by the chemical potential shift on the 
hole (electron)-doped side. This is consistent with our 
observation of different behaviors of the core level shift on 
hole- and electron-doped sides. We note that the change 
of the Madelung term AVm may not be same on hole- 
and electron-doped sides, which can further enhance the 
difference of the core level shift on the two sides. 
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FIG. 4: (a) Phase diagram of the hole- and electron-doped 
Ba-122 systems taken from references ^21^ and ^22^, respec- 
tively. Tc, TsDW and Tstr refer to the superconducting, the 
SDW and the tetragonal to orthorombic structural transi- 
tions, respectively, (b) Doping dependence of the Lindhard 
function at the M point (near-nesting wave vector) normal- 
ized by its value at the zone center. The Lindhard function 
was obtained by using LDA calculation. 



The above analysis suggests that a rigid band picture 
constitutes a good first approximation of the evolution of 
the chemical potential in the 122 family of iron-pnictides. 
We now ask a simple but fundamental question: is FS 
near-nesting able to explain the electron-hole asymme- 
try of the superconducting domes shown in the phase 
diagram of the 122-pnictides of Fig. 4(a)? To answer 
this question, it is necessary to compute the spin suscep- 
tibility. It is especially important to understand how the 
susceptibility evolves at the near-nesting (or AF) wave 



vector. We use the band structure calculated by LDA 
to compute the doping dependence of the Lindhard spin 
susceptibility at the near- nesting wave vector [23l [24] . 
We limit our calculations to the elastic component of the 
spin susceptibility. The results are displayed in Fig. 4(b). 
Interestingly, the hole- and electron-doped sides exhibit 
a strong asymmetry: while the Lindhard function de- 
creases monotonically on the electron-doped side (with a 
small shoulder around ~0.12), it keeps a high value for a 
wide hole doping range before starting to decrease. Re- 
markably, the maximum value of the calculated suscep- 
tibility is obtained near the experimental optimal hole 
doping, and the Lindhard function using the FS nest- 
ing wave vector tracks the superconducting transition 
qualitatively well. It is important to note that such an 
asymmetry in the Lindhard susceptibility would lead to 
a higher dielectric function and consequently to a larger 
screening effect on the hole-doped side. This is qualita- 
tively consistent with the observed larger difference be- 
tween the core level shift and the chemical potential shift 
on the hole-doped side shown in Fig. 3(a). We caution 
that the non-magnetic LDA calculations are no longer 
valid in the spin density wave (SDW) state because the 
band structure undergoes unconventional band folding 
that leads to the formation of Dirac cones [25] . 

The basic reason for electron-hole asymmetry in the 
calculated Lindhard function is as follows. The effec- 
tive masses of the holelike bands, especially the /3 band, 
are larger than that of the electronlike bands at the M 
point, as observed by ARPES [18 and quantum oscil- 
lation experiments [26]. To satisfy the Luttinger theo- 
rem, their top of band at zero doping must thus be closer 
to Ep than the bottom of the electron bands. Indeed, 
even for optimally hole-doped samples, the top of the 
a band is located only 25 meV above Ep [TO]. As a 



consequence, the holelike bands sink below Ep with elec- 
tron doping much faster than the bottom of the electron 
bands are pushed above Ep with hole doping. There- 
fore, the FS near-nesting conditions are more robust in 
the hole-doped case. The built-in asymmetry regarding 
the FS near-nesting condition on the electron- and hole- 
doped sides offer a simple but powerful clue that the FS 
near-nesting with the AF wave vector triggers supercon- 
ductivity in the pnictides. 

In conclusion, we have presented the doping depen- 
dence of the chemical potential in the 122 family of iron- 
pnictides. As a first approximation, our results are con- 
sistent with a rigid band shift and with renormalized 
LDA calculations. The doping dependence of the As 3d 
core levels does not follow that of the chemical poten- 
tial, suggesting a non- negligible screening effect. Within 
the rigid band shift approximation, the calculated Lind- 
hard function at the FS-nesting wave vector based on the 
LDA band structure reveals an electron-hole asymmetry 
in the iron pnictides, which matches well with the ob- 
served electron-hole asymmetry of the superconducting 
domes in the phase diagram. Our findings reveal the im- 
portance of FS near-nesting in the pairing mechanism of 
the iron-based superconductors. 
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